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1. Introduction 


Following the realisation that a simple iterative strategy for bringing the 
flexible walls of two-dimensional test sections to streamline contours was too slow for 
practical use 1 , Judd proposed 2 ^ developed and placed in service 4 what was, as far as we 
know, the first Predictive Strategy (sometimes called a "one-step” method, but see 
comments on the use of this phrase in Section 3). This Strategy, built into a tunnel’s 
control system, makes use of measurements at the flexible walls of the test section to 
predict the magnitudes of the adjustments to their shapes required to eliminate their 
interferences at the model. 

During the following years (1976 to date) the software was further developed 5 * 6 
and extensively used and proved up to transonic speeds 7 * 11 . The later developments 
described in Section 4 were in the form of refinements and did not involve any change in 
the underlying principles. 

The Predictive Strategy reduced by 75% or more the number of iterations of wall 
shapes, and therefore the tunnel run-time overhead attributable to the streamlining 
process, required to reach satisfactory streamlines. As a matter of policy the Strategy has 
been used to eliminate, as far as is experimentally possible the top and bottom wall 
interferences. However it should be noted that as a means for reducing the streamlining 
run-time overhead there remains the option of compromise in the quality of the 
streamlining coupled with the application of modest corrections. 

Because the Strategy is rapid and well proven, it is felt that it would be useful to 
give a detailed description of the software for others easily to adopt. The Strategy works 
well in two-dimensional testing at any set of conditions up to those which result in the 
airfoil’s shock just extending to a streamlined wall (usually this would be the suction 
surface shock just extending to the nearest wall) in a suitably designed test section 1 8 

The Strategy was first implemented in software associated with the running of 
the low speed Self Streamlining Wind Tunnel (SSWT) in 1976, then on the fully 
automated Transonic Self Streamlining Wind Tunnel (TSWT) in 1979, both at the 
University of Southampton, U.K., where the software is still available for use in routine 
two-dimensional testing 910 *! 1 . The simplifications and approximations in its theoretical 
formulation were influenced by the limited computing power available to the team in the 
early days, and the algorithm happened also to have been programmed first in BASIC, 
both influences still being visible in the software. More recently the software has been 
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installed in a computer which controls the flexible walls of the adaptive walled test section 
in the 0.3-m Transonic Cryogenic Tunnel at NASA Langley Research Center. 

The Strategy utilises the velocity distributions along both sides of each flexible 
wall sketched in Figure 1. The real-side velocity distributions are calculated from 
measurements of static pressures along the insides of the walls, while the velocities on the 
outsides of the walls, generated by the imaginary flowfields, are derived by calculation 
using data from the preceding run. This preceding run may have been the preceding 
iteration in a series performed with a particular model, but in fact the wall shapes and 
corresponding imaginary side velocity distributions derived from any previous run may be 
set and used* . The Strategy makes use of this wall information in predicting new wall 
contours which will eliminate the combined top and bottom wall interference present 
during the current run, while simultaneously providing the imaginary-side velocity 
distributions over the new contours 

Wall loading is the evidence of interference: if the real and imaginary velocities 
differ at any point along a wall then the wall is loaded at that point (and in general this is 
so everywhere) and therefore the line followed by the wall is not that of a streamline in the 
infinite fiowfield. The object of the Strategy is to predict the wall movement required to 
eliminate the loading and therefore the interference. Then the wall will be streamlined 

The procedures of the Strategy are embodied in the FORTRAN subroutine WAS 
(standing for Wall Adjustment Strategy) which is written in a general form. The 
following sections of the report begin with a brief description of the essentials of the test 
section hardware, followed by the underlying aerodynamic theory which forms the basis of 
the Strategy. The subroutine is then presented as the Appendix, broken down into 
segments with descriptions of the numerical operations underway in each, with 
definitions of variables. 

Two points should be noted. Firstly, the flexible walls need to be adjusted for 
constant Mach number when the test section is empty to allow for the growth of boundary 
layers, giving what is called ~aerodynamically straight ” wall shapes. The shapes are 
functions of Reynolds and Mach numbers, and are not set with the aid of WAS because in 
moving the walls in the desired direction the subroutine introduces perturbations into the 


f The word * run ” is used here in the context of data gathering: a run is a period during 
which wall pressures (and perhaps other data) are being gathered. 
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imaginary flowfields which should not exist at this stage in the use of the tunnel. 
Secondly, when streamlining around a model, the wall contours which are set must allow 
for the variations in the displacement thicknesses of their boundary layers which are 
induced by the model's pressure field. 


2. Essentials of test section hardware 

The test section comprises a pair of rigid sidewalls which support the model in 
two-dimensional testing, and top and bottom walls made from a convenient flexible 
material. The flexible walls are fitted with a number of jacks which allow the shapes to be 
controlled. The walls are bent by the jacks in single curvature only, are cantilevered at 
their upstream ends and, for minimum interference from length-truncation effects 3 , are 
relatively long and symmetrically disposed fore and aft of the centre of lift of the model. 
The spacing of the jacks need not be regular; in fact it is usual practice to pitch the jacks 
more closely in the region of the model than elsewhere because of the stronger curvature 
in that region. 

The wall streamlining process described here relies on measurements of the 
positions of the walls at each jacking point together with measurements of the wall 
centreline static pressures, also at each jacking point. Reference Mach number is derived 
from reference pressures measured in the usual way at the upstream end of the test 
section. 


3. Basic theory of the strategy 

In its basic form Judd's Predictive Wall Adjustment Strategy applies to the case 
of a single impervious thin wall and a model, both lying in an otherwise undisturbed 
infinite flowfield. His theory applies to the general case of the unstreamlined wall the 
shape of which is known together with the velocity distributions along each side. There is 
no assumption of prior knowledge of the aerodynamic behaviour of the model, neither are 
model measurements a necessary adjunct to the streamlining process. 

The wall is loaded as it does not yet follow the desired line of an unloaded 
streamline in the infinite flowfield. Manifestations of the loading are the differing 
pressures and associated velocities on either side, although the latter are used for 
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convenience in this section. The physical presence of the wall and the distribution of 
velocity difference across it may be replaced by a notional vorticity distribution at the 
wall. The velocity jump (between that on the real side and that on the imaginary side) is a 
direct measure of the local strength of the vorticity. The distribution of vorticity has the 
characteristic that the velocity component induced by it in a direction normal to the wall 
just cancels the sum of the components from other sources thus preventing through-flow. 
One other source of normal velocity component is the model. 

The situation which has just been described, that is a requirement for the 
distribution of vorticity to prevent through-flow at a point, is eliminated by changing the 
slope of the wall at the point so that the vorticity *s contribution to through-flow is replaced 
by a change in the component from the free stream. This is done along the whole wall to 
remove the vorticity everywhere. The operation will be perfect provided the other sources 
of normal velocity component remain unchanged. 

At streamwise station x along a wall positioned above the model the difference in 
velocity is represented by local vorticity of strength 


r( X ) = u(x) - vu) 


where U(x) is the real-side velocity distribution (derived from pressure measurements) 
and V(x) is the imaginary-side (calculated) velocity distribution. A velocity component 
v( 0 normal to the wall at streamwise station £ is induced by the distribution of vorticity. 
For small slopes this is approximately given by the integral of the elemental contributions 
of vorticity atx: 


v(f) = 


2n 


[" n.t )dx 

I (£-x) 


The slope of the wall is adjusted by the amount which is required for a change in the 
normal component of the free stream velocity to just oppose that due to the vorticity. This 
requires (for small values of slope) an increment in slope which is given by the 
approximation 
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or 
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r (x)<f* 

(x-0 


A 




2nU 

X) 


where f/® = free stream velocity. This in turn is integrated leading to the change in wall 
deflection Ay((). 

Following the removal of the vorticity there are adjustments to velocity either 
side of the wall amounting to half of the imbalance existing before movement. This is the 
method by which the velocity distribution is derived for the imaginary side of the wall 
shape which is to be set for the next run. The increment in imaginary-side velocity at 
station x arising from the elimination of the vorticity amounts to (U(x) - V(x))/2. Hence 
the imaginary side velocity for the new shape of wall is 


V(x) + 


U(x) - V(x) 


This basic theory appears to offer immediate streamlining, the so-called one-step 
method. However, a one-step method which does not invoke a knowledge of the model’s 
aerodynamic behaviour would require the behaviour not to change with wall shape, 
whereas the whole of adaptive-wall work arises because model behaviour is dependent on 
test-section boundary conditions, in this case wall shape. A further change to the test 
section flow arises from the second wall which is being streamlined simultaneously 
These interaction mechanisms cause the wall’s predicted shape not to correspond to the 
required streamline. The modifications to the strategy to account for these effects are 
introduced in the next section. 


4* Modifications in service 
4.1 Coupling and scaling 

If there were no other changes in the flowfield then the wail loading could be 
expected to become zero with the new shape. However there will be a change in the 
behaviour of the model induced by the movement of the wall, but more importantly the 
requirement to adjust the opposite wall to bring it also to zero loading will introduce a 
strong interaction. The simultaneous adjustment of each wall using the above simple 
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algorithm does not lead to convergence of the walls to streamlines. Allowance must T>e 
made for what may be regarded as, for long wavelength components of wall movement, a 
one-dimensional continuity effect, a strong source of coupling. Convergence can be 
obtained by feeding a proportion of the demanded movement of one wail to the other. The 
process is now iterative because of the wall-model-wall interactions and, in this form, the 
software results also in an overshoot. That is, the predictions of wall movement are 
somewhat exaggerated. The latter is reduced by scaling down the predicted wall 
movements before accounting for the coupling effect. Empirically determined coupling 
and scaling factors are used. For each of these modifications to wall shape there are 
appropriate adjustments to the calculation of imaginary-side velocities. 

4.1.1 Scaling 

Factoring the distribution of vorticity along a wall by factor SF results in the 
same factoring of slope, wall deflection and increment in imaginary-side velocity. The 
scaled imaginary-side velocity V s at station x along the next wall contour to be set is then 
given by 


V ( X ) = V(x) + SF. 

s 


U(x) ~ V(x) 
2 


for the top wall and similarly for the bottom. 

4.1.2 Coupling 

Coupling requires a proportion CPLF of one wall’s movement to be implanted in 
the other. The modification to the shape (and slope) of the wall receiving the implant 
introduces increments to its imaginary-side velocity distribution. The increment is 
identical to the increment in velocity on the opposite wall had the opposite wall itself been 
moved an amount factored by the coupling factor, in just the same way as when scaling 
Hence the coupled imaginary side velocity V c (x) for one wall is given by 


V (x) 


= V (x) 

$ 




+ CPLF V (x) - L r (x) 
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where suffix V denotes velocities over the opposite wall. The adjustments to the 
imaginary velocities arising from scaling and coupling are carried out simultaneously for 
both walls. Typical values of the factors for both walls are: 

coupling 0.35 

scaling 0.8 

4.2 Checking of imaginary-side velocities 

A key issue is the accuracy of the imaginary-side velocity predictions, since the 
choice of wall shapes and the judgment of whether or not they are streamlined depend on 
the predictions. This issue has been addressed in several ways. Firstly the validity (in 
terms of introducing errors of acceptably small size) of certain approximations in the 
theoretical basis of the Strategy has been investigated 3 with the conclusion that the errors 
are compatible with those arising from other sources, for example experimental error. 
Therefore the reduction of the computational complexity, inherent in the use of a 
simplified algorithm, was justified. Further checks on the velocities predicted by the 
Strategy have included: 

- analytic checks using straight and streamlined wall information derived 
from potential flow theory. 

- the use of source-sink representations of wall shapes to compute the 
imaginary-side velocities < 1.5.12) 


- high subsonic verification of imaginary-side velocity distributions with a 
streamline curvature program. 

- experimental verification of velocity by building a top wall contour into 
the bottom wall of an empty test section then re-streamlining the top 
wall and measuring the bottom wall real-side velocity distribution < b8). 

Each of these has led to the conclusion that the imaginary-side velocities 
computed by the Strategy are reliable. 

4.3 Compressibility 

The wall adjustment strategy of Section 3 is based on potential flow theory, but 
linearised compressible flow corrections were introduced in the following manner to allow 
testing at high subsonic speeds. The various tunnel pressure measurements, in terms of 
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pressure coefficients, C p c are converted to their equivalent incompressible coefficients C p \ 
using 


C = PC where p = V l — M 2 
pi pC r 30 

and Ma> is the reference Mach number. 

Velocities derived from the incompressible pressure coefficients are utilised in 
the strategy. The predicted wall movements and the imaginary-side potential flow 
velocities are stored, available if required for futher iterations. Some earlier 
publications^ contained an error, a factoring of movement demands by P, which is 
corrected here. 

This extension to the Strategy successfully allowed testing up to speeds just 
giving sonic flow at one of the flexible walls. 
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APPENDtX: SOFTWARE DESCRIPTION 


The Predictive Strategy is contained in the subroutine WAS (Wall Adjustment 
Strategy), the essentials of which are presented in the following abstract. There is no 
specific reference to any particular test section and the subroutine may be used in the 
control of any test section as long as it is of the type described in Section 2 and is used in 
two-dimensional testing. 

The following data inputs are required by the subroutine: 
i) Variables: 


NOCPT Number of computing points along one flexible wall (Top 

and bottom walls assumed the same). The computing 
points are located at each of the streamlining jacks and at 
two dummy positions upstream and downstream as 
indicated on Figure 2. The dummy positions are referred 
to as dummy jacks. 


FMACH 

PSTATIC 

TWSF,BWSF 

TWCPLF,BWCPLF 


Free stream Mach Number 
Reference Static Pressure, Pref 
Top, bottom wall scaling factors 
Top, bottom wall coupling factors 


ii) Arrays: 

XJACK(*) Longitudinal co-ordinates of the wall computing points. 

The origin of the wall co-ordinates is the wall anchor point. 

TOPWP(*),BOTWP(*) Static absolute pressures at the wall computing points 

obtained from measurements along the centrelines. 
Dummy wall pressures are assumed: equal to P re f at the 
upstream positions and equal to that at jack N for the 
downstream positions. 
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TWVEL(*),BWVEL(*) 


Normalised velocity perturbations at the computing points 
in the imaginary flows over the current wall contours. 
Normalising is relative to free stream velocity. 


During the current run, the flexible walls have been set to known contours. By 

this we mean that the contours and their imaginary-side velocity distributions are known. 

The data is input from file. 

Data output includes the arrays: 

TWMOV(*),BWMOV(*) Required movements of top, bottom wall jacks from their 

current positions for streamlining. The convention is 
positive upwards for both walls. 

TWNVEL(*),BWNVEL(*} Top, bottom wall imaginary-side normalised velocity 

perturbations at computing points, which will apply to the 
next contours to be set (using TWMOV, BWMOV). 


The software is now broken down into logical segments associated with 
numerical procedures. In this code the word 'velocity 'represents a velocity normalised by 
the free stream velocity. 

Segment 1 

In this segment the real-side wall pressures at computing points, measured 
during the current run, are converted to pressure coefficients, then to equivalent 
incompressible coefficients. The velocity differences between real and imaginary flows 
then lead, after scaling and coupling operations, to the new velocities which will exist on 
the imaginary sides of the next contours to be set. 

Notes on code: 

Line 52 BETA = Prandtl Glauert Factor = V l - M*, 2 

Line 53 Q1 = Dynamic Pressure = y P re f M*r/2 
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Lines 54-84 

Loop calculating velocity data at each wall computing point. 




[I = 1 to NOCPT1 



Line 55 

TCPC 

= Measured top wall pressure coefficient. 



Line 59 

TCPI 

= Equivalent incompressible top wall pressure coefficient. 

Lines 60-61 

TVEL 

= Current real-side velocity perturbation on top wall (“/IT*) 



at computing point I. 



Line 62 

TWVDIF(I) 

= Velocity imbalance across the top wall 



Line 63 

TWVSQ(I) 

= Top wall real-side velocity squared (used in assessment of 



streamlining quality, not used in the strategy). 


Lines 64-72 

Repeat of lines 55-63 for the bottom wall computing point I. 


Lines 76,77 

TWN'VEL(I),B\VN'VEL(I) = Top, bottom wall imaginary-side 

new 


velocity perturbations with scaling. 





= v cn = uin + sf . 

s 

UU] - VT/1 
2 


Line 81 

TNVEL 

= Working store for TWNVEL(I). 



Lines 82,83 

Adjustment of top, bottom wall velocities for coupling 

TWN'VEL(I) and 


BWNVEL(I) now represent the velocity perturbations which will exist on 
the imaginary-sides of the contours to be called by TWMOV( ) and 
BWMOV( ) respectively in Segment 4. 

TWNVEL(I), BWN'VEL(I) = V C [I] 

= V S [I] + CPLF(V 3 [I] o -U[I1 0 ) 
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Subroutine WAS Software Listing 
Segment l 


47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 


C 

C 

C 

C 

C 
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C 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 
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COMPUTE THE VELOCITY IMBALANCE/WALL VORTICITY AT EACH 
WALL COMPUTING POINT AND THE EXTERNAL VELOCITIES 
FOR THE NEXT PREDICTED WALL CONTOURS 

BETA = SQRT ( 1- (FMACH*FMACH) ) 

Q1 = 0.7 * PSTATIC * FMACH * FMACH 

DO 5 I = 1 , NOCPT 

TCPC = (TOPWP (I) -PSTATIC) /Q1 

APPLY PRANDTL-GLAUERT FACTOR TO MEASURED TOP WALL CPS 

TCPI = BETA*TCPC 
TVRATIO = SQRT ( 1-TCPI ) 

TVEL = TVRATIO- 1 
TWVDIF (I) = TVEL-TWVEL ( I ) 

TWVSQ(I) = (TVEL+1) * (TVEL+1) 

BCPC = (BOTWP (I) -PSTATIC) /Q1 

APPLY PRANDTL-GLAUERT FACTOR TO MEASURED BOT. WALL CPS 

BCPI = BETA*BCPC 
BVRATIO = SQRT ( 1-BCPI ) 

BVEL = BVRATIO- 1 

BWVDIF (I) = BWVEL(I) -BVEL 

BWVSQ(I) = (BVEL+1) * (BVEL+1) 

APPLY SCALING FACTORS TO THE EXTERNAL VEL. CALCULATIONS 

TWNVEL ( I ) = TWVEL ( I ) + (TWSF*TWVDIF ( I ) /2 ) 

BWNVEL(I) = BWVEL(I) - (BWSF*BWVDIF ( I ) /2 ) 

APPLY COUPLING FACTORS TO THE EXTERNAL VELOCITIES 

TNVEL = TWNVEL ( I ) 

TWNVEL ( I ) = TWNVEL ( I ) + ( BWCPLF* ( BWNVEL ( I ) -BVEL) ) 

BWNVEL ( I ) = BWNVEL (I) + (TWCPLF* (TNVEL-TVEL) ) 

CONTINUE 
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Segments 2 and 3 


The velocity component normal to a point on a wail induced by its longitudinal 
distribution of vorticity is calculated. This velocity component is then used to determine a 
local change of wall slope. With one exception per wall the points on a wall where the 
induced velocity is determined lie mid-way between computing points, the so-called 
mid-jack points. A cubic is fitted through four adjacent vorticity data points beginning at 
the upstream computing point, and the velocity induced at a mid-jack point by the central 
patch of vorticity (lying between the second and third of the vorticity data points) is 
determined analytically. The group to which the cubic is fitted is then moved one 
computing-point downstream and the process repeated until the last computing point is 
reached, summing the contributions to induced velocity at the same mid-jack point. In 
Segment 2 the streamwise coordinates of the mid-jack points are first determined, then 
the coefficients in the cubic equations for all curve fits. Segment 3 executes the 
integrations and finally determines the required increments in the wall slopes at all 
mid-jack points. 

Notes on code: 


Segment 2 

= The longitudinal co-ordinate of just the first mid-jack 
point is made to coincide with the wall anchor point 
because this point has zero slope and deflection. 

= The longitudinal co-ordinate of the 1th mid-jack point 
lies mid-way between XJACK(I) and XJACK(I + 1). 
NCPT1 = NOCPT-2. 

Lines97-160 Loop for top and bottom wall calculations (NN = 1 and 2 respectively) 

Lines 102-126 Loop to compute coefficients for piecewise cubic curve fit to the wall 
vorticity (in the form r(x >fU^) at each computing point. 
NCPT2 = N'OCPT-3. 

Lines 104-110 Load sets of four XJACK values and BWVDIF or TWVDIF values into 
arrays X and VEL respectively. 


Line 90 XMIDJ(l) 


Lines 94-95 XMIDJ(I) 
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Lines 111-125 


A curve is fitted to the four sets of vorticity data in arrays X and VEL 
and the coefficients of the cubic ax 3 + bx 2 +■ cx + d computed, where 

a = CUBCOE (IL,4) 
b = CUBCOE (IL,3) 
c = CUBCOE (IL, 2) 
d = CUBCOE (IL.l) 


where IL is the counter for each patch of wall vorticity. The geometry of 
the curve-fit is illustrated on Figure 3. 
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Segment 3 
Lines 133-159 

Line 134 

Line 137 

Lines 138-154 

Line 139 
Line 144 
Line 153 

VELSUM = 
Lines 156,158 


Loop to integrate the vorticity along each wall at the mid-jack points 
numbered from 2 to NOCPT-2. 

X0 = Co-ordinate of mid-jack point for which each integration is 
made. 

VELSUM = Total sum of vertical velocity induced by the complete 
wall vorticity. Initially set to zero. 

Loop to perform the analytical integration of the vorticity-induced 
normal velocity at mid-jack point Xo for I = 1 to NOCPT-3, where I is 
the counter for each patch of wall vorticity. See Figure 3. 

XI = Co-ordinate of lower limit of integration = XJACKfl + 1). 

X2 = Co-ordinate of higher limit of integration = XJACKd + 2). 

VELSUM = Summing operation for vertical velocities induced by each 
patch of wall vorticity where 


( x2 ( COEFFti + COEFFl.x + COEFF2.X 2 + COEFFZx x 
y£LS(/Af+ ; dx 

Li V (x — Xu) J 


This is solved analytically using four standard integrals coded in lines 
147 to 152. Since X 0 * xl or x2 the singularity x = X 0 is avoided 

TSLOPE(J),BSLOPE(J) = VELSUM/2n = top, bottom wall required 

change in local wall slope. 
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Subroutine WAS Software Listing 
Segment 2 


86 c 

87 C MAKE THE FIRST MID-JACK CO-ORD AT THE WALL ANCHOR POINT 

88 C TO ENSURE A ZERO WALL SLOPE AT THIS LOCATION 

89 C 

90 XMIDJ(l) = XJACK(2 ) 

91 C 

92 C DETERMINE OTHER MID- JACK CO-ORDS BETWEEN WALL COMPUTING 

93 C POINTS 

94 DO 15 I = 2 , NCPT1 

95 XMIDJ(I) = (XJACK(I)+XJACK(I+1) )/2 

96 15 CONTINUE 

97 DO 25 NN = 1,2 

98 C 

99 C PIECEWISE CUBIC CURVE FIT TO THE WALL VORTICITY USING 

100 C SETS OF FOUR COMPUTING POINTS (LABELLED 1,2, 3 ,4) 

101 C 

102 DO 95 IL = 1 , NCPT2 

103 I = IL - 1 

104 40 DO 35 J = 1,4 

105 X(J) = XJACK ( I+J) 

106 IF (NN. EQ. 1) GO TO 50 

107 VEL(J) = BWVDIF ( I+J) 

108 GO TO 35 

109 50 VEL(J) = TWVDIF ( I+J) 

110 35 CONTINUE 

111 V0 = ( VEL ( 3 ) -VEL (2))/(X(3) -X ( 2 ) ) 

112 VI = VEL ( 2 ) -V0*X ( 2 ) 

113 DIST1 = 1/ (X (4 ) -X ( 1) ) 

114 V2 = (VEL (4 ) -V0*X (4 ) -VI) / ( (X ( 4 ) -X (2) ) * (X ( 3 ) -X ( 4 ) ) ) 

115 V3 = (VEL (1) -V0*X (1) -VI) / ( (X ( 1) -X (2))*(X(3) -X ( 1) ) ) 

116 V4 = DIST1* (V2-V3) 

117 V5 = V3-V4*X (1) 

118 DIST2 = X ( 2 ) + X ( 3 ) 

119 C 

120 C CALCULATE COEFFS. FOR EACH PIECEWISE CUBIC CURVE FIT 

121 C 

122 CUBCOE ( IL, 1 ) = Vl-X ( 2 ) *X ( 3 ) *V5 

123 CUBCOE ( IL, 2 ) = V0+V5*DIST2-V4 *X (2 ) *X ( 3 ) 

124 CUBCOE (IL, 3) = V4*DIST2-V5 

125 CUBCOE (IL, 4) = -V4 

126 95 CONTINUE 


16 



Subroutine WAS Software Listing 
Segment 3 


127 C AT EACH MID-JACK PT ., INTEGRATE THE VORTICITY ALONG EACH 

128 C WALL TO FIND THE INDUCED VERTICAL VELOCITIES, ASSUMED 

129 C NORMAL TO THE TOP AND BOTTOM WALLS, WHICH MUST BE 

130 C CANCELLED BY CHANGES IN THE FREE STREAM COMPONENT 

131 C CAUSED BY LOCAL ADJUSTMENT OF WALL SLOPE 

132 C 

133 DO 45 J = 2 ,NCPT1 

134 X0 = XMIDJ(J) 

135 XOSQ = X0*X0 

136 XOCUB = X0SQ*X0 

137 VELSUM =0.0 

138 DO 55 I = 1 , NCPT2 

139 XI =XJACK( 1+1) 

140 COEFFO = CUBCOE (1,1) 

141 COEFF1 = CUBCOE (I, 2) 

142 COEFF2 = CUBCOE (I, 3) 

143 COEFF3 = CUBCOE (I, 4) 

144 X2 = XJACK(I+2) 

145 X2SQ = X2 * X2 

146 X1SQ = XI * XI 

147 SUMO = COEFFO+COEFF1*XO+COEFF2* (XOSQ) +COEFF3* (XOCUB) 

148 X3 = ABS (X2-X0)/ABS (X1-X0) 

149 X4 = ALOG (X3) 

150 SUM1 = (COEFF1+COEFF2*XO+COEFF3*XOSQ) * (X2-X1) 

151 SUM2 = (COEFF2+COEFF3*XO) * ( (X2SQ) - (X1SQ) ) /2 

152 SUM3 = COEFF3 * ( (X2SQ*X2 ) - (X1SQ*X1) ) /3 

153 VELSUM = VELSUM+SUM0*X4+SUM1+SUM2+SUM3 

154 55 CONTINUE 

155 IF (NN.EQ.2) GO TO 60 

156 TSLOPE(J) = VELSUM/ 6. 28319 

157 GO TO 45 

158 60 BSLOPE(J) = VELSUM/6 . 28319 

159 45 CONTINUE 

160 25 CONTINUE 
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Segment 4 

The increments required in wail slope are available at each mid-jack point. 
These are now integrated to provide wail movement, the integration beginning at the 
anchor point which remains fixed with zero slope. The general technique is to fit the 
quadratic equation ax 2 + bx + c through three adjacent values of wail slope increment (as 
a function of streamwise position). This quadratic equation is then integrated giving a 
cubic which passes through the predicted changes of wall positions of each of the three 
mid-jack points shown on Figure 4. The first three coefficients of the cubic equation 
Ax 3 + Bx 2 + Cx + D are related to those of the quadratic equation as follows 


A = 


a 

3 * 



C = c 


The integration is performed between the x-limits of the two jacks which are straddled by 
this group of mid -jack points, giving the relative change of curve (wall) height between the 
two jacks. The process is repeated step-by-step along the whole test section from the fixed 
upstream end, giving the required movement (and shape) of the complete wail. 

Notes on code 


Lines 164,165 Initialise top, bottom wall jack movement integrands. 

Lines 169,170 TSLOPE(l),BSLOPE(l) Top, bottom wall slopes at XMIDJ(l) (the wall 

anchor points) set to zero. 


Lines 176-216 Loop calculating the jack movement demands for wall streamlining at 
computing point (I + 2), namely Jack I, NCPT3 = NOCPT-4 


Lines 187-195 Determination of cubic coefficients for top wall position changes, where 

a b 

A = - ; B = - ; C = c 
3 2 


Line 196 TMOV = Required movement of jack I on top wall. 

Lines 197-204 Determination of bottom wall coefficients. 
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Line 205 
Lines 209-215 


BVOV = Required movement of jack I on bottom wall. 

Scaling , then coupling of top and bottom wall jack movement demands. 
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Subroutine WAS Software Listing 


161 C 

162 C 

163 C 

164 

165 

166 C 

167 C 

168 C 

169 

170 

171 C 

172 C 

173 C 

174 C 

175 C 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 C 

191 C 

192 C 

193 

194 

195 

196 

197 

198 

199 C 

200 C 

201 C 

202 

203 

204 

205 

206 C 

207 C 

208 C 

209 

210 

211 C 

212 C 

213 C 

214 

215 

216 65 


Segment 4 

INITIALISE WALL MOVEMENT DEMAND ACCUMULATORS 

TMOV =0.0 
BMOV =0.0 

SET WALL SLOPES AT THE WALL ANCHOR POINTS EQUAL TO ZERO 

TSLOPE ( 1) =0.0 
BSLOPE ( 1) =0.0 

FIND THE JACK MOVEMENTS REQUIRED FOR WALL STREAMLINING, 
BY PERFORMING INTEGRATIONS OF PIECEWISE QUADRATIC 
CURVES FITTED TO SETS OF THREE WALL SLOPES 

DO 65 I = 1 , NCPT3 

11 = 1+1 

12 = 1+2 

TSGRAD = (TSLOPE (12 ) -TSLOPE ( II) ) / (XMIDJ ( 12 ) -XMIDJ (II)) 
BSGRAD = ( BSLOPE ( 12 ) -BSLOPE (II) ) / (XMIDJ ( 12 ) -XMIDJ (II)) 
XJ1SQ = XJACK(Il) *XJACK(I1) 

XJ2SQ = XJACK (12 ) *XJACK ( 12 ) 

XJ1CUB = XJ1SQ * XJACK (II) 

XJ2CUB = XJ2SQ * XJACK (12) 

XI = XMIDJ (I) -XMIDJ (II) 

X2 = XMIDJ (12) -XMIDJ (I) 

PI = (TSGRAD- (TSLOPE (I) -TSLOPE (II) ) /XI) /X2 
P2 = TSGRAD - PI* XMIDJ ( 12 ) 

X3 = XJACK (12) -XJACK (II) 

TOP WALL - MOVEMENT DEMAND CUBIC COEFFICIENTS 
A = PI/ 3 

B = (P2-P1*XMIDJ(I1) )/2 
C = TSLOPE (II) - P2*XMIDJ (II) 

TMOV=TMOV+ (A* (XJ2CUB-XJ1CUB) ) + (B* (XJ2SQ-XJ1SQ) ) + (C*X3 ) 
PI = (BSGRAD- (BSLOPE (I) -BSLOPE (II) ) /XI) /X2 
P2 = BSGRAD - PI* XMIDJ ( 12 ) 

BOTTOM WALL - MOVEMENT DEMAND CUBIC COEFFICIENTS 
A = PI/ 3 

B = (P2-P1*XMIDJ(I1) )/2 
C = BSLOPE ( II ) - P2*XMIDJ(I1) 

BMOV=BMOV+ (A* (XJ2CUB-XJ1CUB) ) + (B* (XJ2SQ-XJ1SQ) ) +(C*X3) 

SCALE JACK MOVEMENT DEMANDS 

STMOV = TWSF * TMOV 
SBMOV = BWSF * BMOV 

COUPLE JACK MOVEMENT DEMANDS 

TWMOV(I) = STMOV+ (BWCPLF*SBMOV) 

BWMOV(I) = SBMOV+ (TWCPLF*STMOV) 

CONTINUE 
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FIC. 1 DIVISION OF THE INFINITE TWO-DIMENSIONAL FLOWFIELD INTO REAL AND 
IMAGINARY PARTS. 


Pressures assumed Pressures assumed 



I Jack + wall static pressure tapping 

V Dummy jack . 

A Anchor point. Slope zero. 

Dummy straight wall extensions 


FIC. 2 REPRESENTATION OF A FLEXIBLE WALL IN THE WALL ADJUSTMENT 
STRATEGY. 
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77 = ax 3 + bx 2 + cx + d 
u oo 



Mid-jack 

point 


Normalised vertical velocity induced at X0 by the patch 
of vorticity extending from XI to X2 is 

1 f X2 (ax 3 ♦ bx 2 + cx + d) 

2-tt J (x - X0 ) dX 

x=X1 


FIG. 3 PIECEWISE ANALYTICAL TECHNIQUE USED TO INTEGRATE THE 
VORTICITY-INDUCED UPWASH AT MID-JACK POINT X0. 
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Demands for 
wall movement , 



O Movements derived from integration of increments 
in wall slope. 

X Interpolated demands for jack movement. 


FIG. 4 INTERPOLATION OF THE MOVEMENT DEMANDED OF JACK I 
RELATIVE TO ADJACENT UPSTREAM JACK 1-1 
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